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f *^ , Absorbing-boundary-condition method and its applications to nuclear responses 

and breakup reactions are reported. The method facilitates calculations of the 
continuum states in the coordinate space of many degrees of freedom. Properties 
of nuclei near drip lines are discussed. 



1. Absorbing boundary condition (ABC) 

Advances in radioactive beams provide us with a good opportunity to study 
physics of weakly-bound finite-quantum systems. Since excitation spectra 



o 
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above the particle-emission threshold are continuum spectra, theoretical 
analysis requires continuum wave functions. We have recently investigated 
an efficient and comprehensive method of treating the continuum 1 ' 2 ' 3,4 ' 5,6 . 
This is practically identical to the one called "Absorbing Boundary Condi- 
tion (ABC) method" in the chemical reaction studies 7 . The method allows 
us to calculate the continuum wave functions with the outgoing asymptotic 
behavior of many-body systems in a finite spatial region where particles 
interact with each other. 

The essential trick for the treatment of the continuum in the ABC 
method is to allow the infinitesimal imaginary part in the Green's func- 
tion, ie, to be a function of coordinate and finite, ie(r). The e(r) should 
be zero in the interacting region and be positive outside the physically rel- 
evant region of space. Wave functions obtained using the ABC method are 
meaningful only in the interacting region. However, this is enough to solve 
scattering problems of particles which interact by a force with finite range. 
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In order to understand how the ABC method is functioning in later 
applications, it is useful to consider the potential scattering of a particle. 
The scattering wave function is given by 

k +) > = |k> + l,. V\k) = |k> + 1^)- (1) 
hi — tl + ie 

The scattering amplitude, /(O), is usually defined by its asymptotic behav- 
ior 

^ +) (r) - exp(ik • r) + /(fi)^*!^ , ( r _ «,), ( 2 ) 
but can be written in a form 

m) = -^( k 'l^l^k +) ) = -7^2 / drexp(-ik' • r)y(r)^ +) (r), (3) 



where f2 is the direction of k' and |k'| = |k|. Equation (3) implies that 

4 + ' 



the f(Q) can be determined by the scattering wave function, i/;^ + '(r), in 



the interacting region where V =/= 0. In other words, behavior of ipu~^ ( r ) 
outside the interacting region is irrelevant to determination of the scattering 
properties. 

Figure 1 shows calculated phase shifts for a nucleon scattered by a 
square-well potential: 

V(r) = L V ° T r<r °' (4) 

w [ for r > r , 

where Vo = —20 MeV and ro = 3 fm. The upper part of Fig. 1 shows the 
phase shifts calculated with a standard numerical procedure. Namely, we 
solve the radial Schrodinger equations to obtain the regular wave functions 
of s-, p-, and d-waves, and determine the phase shifts at a certain point of 
r > r . 

On the other hand, Eqs. (1) and (3) supply another way of solving the 
problem. The iV'scat) m -^q. (1) satisfies the Schrodinger equation with a 
source term 

(E-H + i e)\4^) = V\k). (5) 

This equation must be solved with the outgoing boundary condition. Here, 
the ABC plays an essential role. The infinitesimal quantity ie in Eq. (5) is 
now replaced by ie(r), then, the equation can be solved with the vanishing 
boundary condition. The conditions and limitations on e(r) are discussed 
in a number of works 1 ' 7,8 . As is emphasized before, the obtained wave 
function, '0iit( r )> is correct only in the interacting region. The scattering 
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Figure 1. Calculated phase shifts for a nucleon scattered by a square- well potential. 
See text for explanation. 

amplitude, f(Q), is calculated using Eq. (3), and arg(/(0)) corresponds to 
the phase shift, which is displayed in the lower part of Fig. 1. The results 
are identical to those of the standard method. The agreement in results 
of the two methods clearly indicates that the continuum in the interacting 
region (r < ro) is properly taken into account in the ABC calculation. 



2. Application to nuclear breakup reactions of Be 

The solution of the two-body problem in Sec. 1 was trivial. The ABC did 
not show any advantage over the standard techniques. However, quantum 
three-body scattering problems are much more difficult, and one can find 
usefulness of the ABC. 

Let us consider a reaction of a projectile, composed of core (C) plus 
neutron (n), on a target nucleus (T). Denoting the projectile-target rela- 
tive coordinates by R and the neutron-core relative coordinates by r, the 
Hamiltonian of this three-body system is expressed as 

H = ~ I^V* + V nC (r) + V nT {Y nT ) + V C t{^ct) (6) 

where fi and m are the reduced masses of projectile-target relative motion 
and neutron-core relative motion, respectively. V n c V n x, Vqt are the 
interaction potentials of constituent particles. 

The wave function may be expressed as a sum of the Coulomb wave in 
the incident channel and the scattered wave. 

*(+>(R,r) = V (+) (R)^o(r) + *scat(R,r) (7) 

where (f>o(r) is the ground state of the projectile, described as a n-C bound 
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state. The \I/ sca t satisfies the following inhomogeneous equation in the ABC, 
{E + e + ie nC {r) + ie PT {R) - H} * scat (R, r) 
= {V nT (v„T) + Verier) - V c }iP (+ Hn)Mr), (8) 

where Vc is the Coulomb distorting potential and eo is the ground- 
state energy of the projectile. One should note that the right hand side 
{V n T + Vct — Vc} ip( + \H)(j>o(r) is a localized function in space, being 
analogous to the fact that right hand side of Eq. (5) is localized. Numerical 
details can be found in our recent paper 6 . We have studied a deuteron 
breakup reaction and compared our results with those of the continuum 
discretized coupled channel (CDCC) calculation 6 . The results agree with 
those of CDCC. In this paper, we report the application to a breakup re- 
action of 11 Be. 

The 10 Be-n potential is taken as Woods-Saxon shape whose depth is 
set so as to produce the 2s orbital binding energy. We adopt the optical 
potential 5 for 10 Be- 12 C and the Becchetti-Greenlees potential 9 for n— 12 C. 
The radial region up to 30 fm and 50 fm are used for R and r, respectively. 
The ie n c( r ) and iepr{R) is non-zero in the region 20 fm < R < 30 fm and 
25 fm < r < 50 fm. The n- 10 Be relative angular momenta are included up 
to I = 3. 

In Fig. 2, we show the elastic breakup cross sections of 11 Be- 12 C reac- 
tion. The filled circles are the result of the ABC calculation and the open 
circles for the eikonal calculation. The elastic breakup cross section is sub- 
stantially larger than that in the eikonal approximation at lower incident 
energy. The failure of the eikonal approximation is apparent at the inci- 
dent energy below 50 MeV/A. There, the quantum-mechanical treatment 
is required for the three-body continuum. 

3. Giant dipole resonances in superdeformed 14 Be 

In studies of giant resonances, effects of the continuum has been treated 
in the random-phase approximation (RPA) with Green's function in the 
coordinate space 10 . However, it is very difficult to directly apply the method 
to deformed nuclei because construction of the Green's function becomes 
a difficult task for the multi-dimensional space 1 . We have shown that the 
ABC method is very useful to treat the electronic continuum in deformed 
systems, such as molecules and clusters 1 . We have also investigated the 
applicability of the ABC in studies of nuclear response calculations 2 ' 3,4 . In 
this section, we discuss a giant dipole resonance (GDR) in 14 Be. 
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Figure 2. The elastic breakup 
cross section in 11 Be- 12 C re- 
action. Quantum calculation 
with ABC (closed circles) is 
compared with the eikonal cal- 
culation (open). 



Figure 3. Calculated El strength in 14 Be. Thin 
solid and dashed lines indicate the response to 
dipole fields parallel and perpendicular to the 
symmetry axis, respectively. Thick line shows 
the total strength. 



We use the ABC in the time-dependent Hartrec-Fock (TDHF) calcu- 
lations on a three-dimensional (3D) coordinate grid. In the real-time cal- 
culations, the linear response is computed by applying the isovector-dipole 
field to the Hartree-Fock (HF) ground state of 14 Be, 

V ext (t)=fcrjl(l-T 2 )e-^}<5(i), (9) 

where k should be small enough to validate the linear response approxima- 
tion. We calculate the expectation values of the El operator as a function 
of time, and then Fourier transforming to get the energy response 1 ' 2 ' 3 ' 4 . 
Since all frequencies are contained in the initial perturbation, the entire 
energy response can be calculated with a single time evolution. 

We use the Skyrme energy functional of EV8 with the SIII parameter 
set 11 . For the time evolution of the TDHF state, we follow the standard 
prescription 12 . The model space is a sphere whose radius is 22 fm. The 
ie(r) is zero in a region of r < 10 fm, while it is non-zero at r > 10 fm. The 
TDHF single-particle wave functions are discretized on a rectangular mesh 
in a 3D real space. Now, we can perform the TDHF simulation together 
with the vanishing boundary condition at r = 22 fm. Time evolution is 
carried out up to T = 10 h/MeV. 

The density distribution of the 14 Be ground state is calculated to have a 
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prolate superdeformed shape. Therefore, we expect the deformation split- 
ting of the GDR peak. The El oscillator strengths of 14 Be are shown in 
Fig. 3. Here, we use a smoothing parameter of T = 1 MeV. The calcu- 
lation predicts the large deformation splitting of about 10 MeV. However, 
the width of each peak is so large that the double-peak structure is almost 
smeared out in the total strength (thick line). This is different from our 
results of stable nuclei 2,4 . We have observed prominent two-peak structure 
in those nuclei. The significant damping width may be a peculiar nature of 
drip-line nuclei. 

4. Summary 

The idea of absorbing boundary condition is presented. The continuum is 
properly treated in the interacting region, although the equations can be 
solved with the vanishing boundary condition. Since we do not need to 
construct the outgoing boundary wave functions, the ABC greatly facili- 
tates the scattering problems for cases of the continuum in many degrees 
of freedom. We apply the ABC to the quantum reaction studies of elastic 
breakup of 11 Be and the linear response calculation for the isovector GDR 
of superdeformed 14 Be. 

This work is supported by the Grand-in-Aid for Scientific Research (No. 
14740146) from the Japan Society for the Promotion of Science. 
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